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Self-gravitating SU(2) Higgs magnetic monopoles exist up to a critical value of the ratio of 
the vector meson mass to the Planck mass, which depends on the Higgs boson mass. At the 
critical value a critical solution with a degenerate horizon is reached. As pointed out by Lue and 
Weinberg, there are two types of critical solutions, with a transition at an intermediate Higgs 
boson mass. Here we investigate this transition for black holes, and reconsider it for the case of 
gravitating monopoles. 

I. INTRODUCTION 

In the SU(2)-Einstein-Yang-Mills-Higgs model, gravitating magnetic monopoles and non-abelian black holes 
exist in a certain region of the parameter space |Q-p[- For a fixed value of the Higgs boson mass, gravitating 
monopoles exist up to a critical value a cr (of the parameter a, which is proportional to the ratio of the vector meson 
mass to the Planck mass), and non-abelian black holes exist up to a critical value a cr (xh) (for small enough values 
of the horizon radius xu |3|). At the critical value a cr a critical solution with a degenerate horizon is reached. In 
particular, for small values of the Higgs boson mass, the critical solution where a horizon first appears corresponds 
to an extremal Reissner-Nordstr0m (RN) solution outside the horizon while it is non-singular inside. 

Recently Lue and Weinberg M reconsidered self-gravitating magnetic monopoles. In particular, they observed, 
that for larger values of the Higgs boson mass, the critical solution is an extremal black hole with non-abelian 
hair and a mass less than the extremal RN value. Exploring the transition between the two regimes, occurring at 
some intermediate value of the Higgs boson mass, Lue and Weinberg were left with a discrepancy between their 
■ analytical and numerical results (J]. 

Here we extend the investigation of Lue and Weinberg to non-abelian black holes, showing the transition to 
persist in the presence of fixed finite event horizon radii Xh- Furthermore, we reconsider the case of gravitating 
monopoles and argue, that the discrepancy between Lue and Weinberg's analytical and numerical results Q can 
be traced back to their limited numerical analysis. 

II. EQUATIONS OF MOTION 

We consider the SU(2) Einstein- Yang-Mills-Higgs action, with a Higgs triplet 
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e is the gauge coupling constant, A is the Higgs coupling constant and v is the Higgs field vacuum expectation 
value. Variation of the action eq. ([!]) with respect to the metric g M „, the gauge field A® and the Higgs field (jf 
leads to the Einstein equations and the matter field equations. 

To construct static spherically symmetric globally regular and black hole solutions we employ Schwarzschild-like 
coordinates and adopt the spherically symmetric metric [^|,^| 

ds 2 = g^dx^dx" = -A 2 Ndt 2 + N^dr 2 + r 2 (dQ 2 + sin 2 9d(/) 2 ) , (4) 

with the metric functions A(r) and N(r), 

, , 2m(r) , . 

N(r) = 1 Li , (5) 

r 

where m(r) ist the mass function. 

For the gauge and Higgs fields we employ the standard spherically symmetric ansatz with vanishing time 
component of the gauge field Jl]-|| 

1 - Kir) - 1 - K(r) 
A t = , A r = , A g = -e — K -L , A^ = e e — — sin6» , (6) 

and 

$=e r H(r)v, (7) 

with unit vectors e r , eg and e^. 

We now introduce the dimensionless coordinate x and the dimensionless mass function fj,, 

x = evr , [i = evm , (8) 

as well as the coupling constants a and /3, 

a 2 = AirGv 2 , (3 2 = 4r . (9) 
9 2 

The tt and rr components of the Einstein equations then yield the equations for the metric functions, 

// = a 2 ( NK' 2 + -Nx 2 H' 2 + ^ + H 2 K 2 + ^x 2 (H 2 - l) 2 ] , (10) 

\ 2 2x 4 / 

and 



A> = a 2 x(^^+H /2 jA, (11) 

where the prime indicates the derivative with respect to x. For the matter functions we obtain the equations 

(ANK'Y = AK(^- T ±+H 2 ^j , (12) 

and 

(x 2 ANH'Y = AH (2K 2 + p 2 x 2 (H 2 - 1)) . (13) 

The equations of motion depend only on the two physical parameters a and j3, eq. (|^). They are related to the 
parameters a and b employed by Lue and Weinberg Q via 

a = 2a 2 = 8nGv 2 , b = U 2 = ( ^-V . (14) 

2 \2mw J 

We note that the embedded RN solution with mass /ioo and unit magnetic charge 

a 2 

fi{x) = Hoo - , A(x) = 1 , K{x) = , H(x) = 1 . (15) 
2x 

is a special solution of these equations. The extremal RN solution, in particular, has horizon Xh, 

x h = x Q = Hoo = a . (16) 
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III. GRAVITATING MONOPOLES 



Let us first consider the globally regular particle-like solutions of the SU(2) EYMH system, the gravitating 
monopoles. Requiring asymptotically flat solutions implies that the metric functions A and [i both approach a 
constant at infinity. We here adopt 

A(oo) = 1 , (17) 

and /i(oo) = /ioo represents the dimensionless mass of the solutions. The matter functions also approach constants 
asymptotically, 

if(oc)=0, iT(oo) = l. (18) 
Regularity of the solutions at the origin requires 

/i(0) = , if (0) = 1 , H(0) = . (19) 

Let us now recall how the gravitating magnetic monopole solutions approach critical solutions, when the vector 
boson mass, i.e. a, is varied, while the ratio of the Higgs boson mass to the vector boson mass, i.e. b, is kept fixed. 
We focus on intermediate values of b. 

Lue and Weinberg 0] realized that there are two regimes of b, each with its own type of critical solution. In 
the first regime b is small. Here the metric function N(x) of the monopole solutions always possesses a single 
minimum. With increasing a the minimum of the function N(x) decreases. For a —> a cr , the critical solution is 
approached, where the minimum of the function N(x) reaches zero at x = xq, eq. jl^). For x > xq the critical 
solution corresponds to an extremal RN black hole solution with horizon radius Xh = xq and unit magnetic charge, 
eq. (15). Consequently, the mass of the limiting solution coincides with the mass of this extremal RN black hole, 
cq. (16). However, in the interior region, x < xq, the critical solution is not singular, in particular N(0) — 1. We 
refer to this approach to the critical solution as RN-type behaviour. RN-type behaviour was seen for gravitating 
monopoles and black holes as well as for gravitating dyons and dyonic black holes |^|| . 

In the newly found B second regime b is large, and the metric function N(x) of the monopole solutions develops 
a second minimum as the critical solution is approached. This second minimum arises interior to the location 
of the first minimum, the solution therefore now possesses an inner and an outer minimum. Once present, the 
inner minimum decreases faster with increasing a than the outer minimum. Consequently, the inner minimum 
reaches zero, when the outer minimum still has a finite value, which is indeed not too different from the value 
it had, when the inner minimum first appeared. Thus for a — > a cr , the inner minimum of N(x) reaches zero at 
x — x* < xq. This second approach to the critical solution is demonstrated in Fig. 1, where we show the function 
N(x) for b = 36.125 and several values of a close to a cr . Here the critical solution possesses an extremal horizon 
at a;* < xq, corresponding to an extremal black hole with non-abelian hair and a mass less than the extremal RN 
value. Consequently we refer to this approach to the critical solution as NA-type (non-abelian-type) behaviour. 

In their paper 0J Lue and Weinberg investigated analytically the transition from RN-type behaviour to NA- 
type behaviour. One of their central analytical results is, that RN-type behaviour is possible only for a > 1.5, 
making a tr = 1.5 the critical value, where the transition from RN-type to NA-type behaviour should occur. This 
lower bound is found by an accurate analytical expansion of the solution about the point x = Xq, where N(x) 
has a double zero. The hypotheses used are mild and, very likely, fulfilled by the solutions. 

When investigating the transition from RN-type to NA-type behaviour numerically Lue and Weinberg fl 
observe, that with increasing b the critical value a cr decreases. For small b the critical value a cr is larger than 
the analytical transition value at r , and the behaviour is of RN-type as expected. As b is increased, and a cr 
decreases towards atr = 1.5 the behaviour remains of RN-type. (a tr = 1.5 would correspond to b m 25.6 01.) 
But surprisingly, as a tT — 1.5 is passed, Lue and Weinberg do not immediately observe the transition to NA-type 
behaviour. Instead they continue to see numerically RN-type behaviour until a value a cr — 1.42, corresponding 
to b m 40, which is clearly below their analytical prediction for the transition. 

Convinced by the analytical argument given by Lue and Weinberg Q for the transition value, we argue in the 
following that the numerically observed discrepancy arises only because of their interpretation of their in accuracy 
limited numerical analysis. 
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Performing an independent numerical analysis (see Appendix for details), we observe that N A- type behaviour 
occurs indeed for lower values of b than b = 40 and thus larger values of a than a — 1.42. Numerically we find 
NA-type behaviour for b > 30. This is illustrated in Figs. 2 and 3, where we present the critical value a cr — 1 
and the ratio x*/xq as functions of b, respectively, in the NA-type regime. These figures should be compared to 
Figs. 10 and 8 of jj), respectively. 

A crucial point of the numerical analysis is illustrated in Fig. 4. When for a fixed value of b in the NA-type 
regime the critical solution is approached, the inner minimum appears at a certain minimal value of a. In Fig. 4 
we present the value of the function N(x) at the outer minimum at this certain minimal value of a (where the 
inner minimum appears) as a function of b. Clearly, for 40 > b > 30, the inner minimum of the metric function 
N(x) appears only, when the value of the function N(x) at the outer minimum is already very small. Indeed, with 
decreasing b it is decreasing from a value on the order of 10 -6 to a value on the order of 10~ 9 . Obviously, with 
decreasing b it becomes numerically increasingly difficult, to observe NA-type behaviour at all. As illustrated 
in Fig. 5 for b — 36.125, once the inner minimum is present, the value of the inner minimum of N(x) decreases 
rapidly with increasing a, whereas the value of the outer minimum remains practically unchanged. 

Even though our numerical analysis is not accurate enough to find NA-type behaviour for b < 30, we expect it 
to be present until atr = 1.5 is reached. Linear extrapolation of the curve in Fig. 2 leads to &tr = 25.58, while (a 
less accurate) extrapolation of the curve in Fig. 4 leads to &tr = 26.7 0. Thus our numerical analysis strongly 
supports the analytical result of Lue and Weinberg Q . 

What concerns their numerical results, we fully agree with their explicit calculations as far as presented in 
However, Lue and Weinberg do not consider numerical solutions, where the value of the outer minimum is smaller 
than 10" 6 f§§. Thus they miss NA-type solutions for values of b smaller than about 40, as seen from Fig. 4. 
Indeed, if we extract from our numerical data a modified "critical" value of a by requiring that the value of the 
minimum of N(x) should be 10~ 6 || instead of requiring that it should be zero, we obtain full agreement with 
Fig. 10 of Lue and Weinberg Q. 



IV. BLACK HOLES 



We now turn to the black hole solutions of the SU(2) EYMH system. Imposing again the condition of asymptotic 
flatness, the black hole solutions satisfy the same boundary conditions at infinity as the regular solutions. The 
existence of a regular event horizon at Xh requires 



Xh 
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fJ-(xh) = 

and A{xh) < oo, and the matter functions must satisfy 

'K 2 - 1 



N'K'\ Xh = K 



and 



2 N'H'\ Xh = H (2K 2 + /3 2 x 2 (H 2 - 1)) | a 



(20) 



(21) 



(22) 



Magnetically charged non-abelian black hole solutions exist only in a limited region of the (a,Xh) plane |l] 0], 
as illustrated for instance for b = in Fig. 7 of 0. When a is varied, while b and the horizon radius Xh are kept 
fixed, magnetically charged non-abelian black holes also approach critical solutions 

Extending the above analysis to the case of black holes, we observe again two regimes of b, each with its own 
type of critical solution corresponding to NA-type and RN-type behaviour, respectively. Since the transition value 
a tr = 1.5 is obtained by performing a Taylor expansion about a point "far" from the horizon and by making use 
of the first few terms of this expansion, we expect this result to hold also in the case of black holes. 

Our numerical results for black holes with horizon radius x^ = 0.01 and 0.1 are included in Figs. 2-4, along 
with the corresponding results for the regular solutions. We observe a complete analogy of the results for regular 
solutions and black holes. With increasing horizon radius Xh, the transition occurs for decreasing values of b. 
Extrapolating the curves of Fig. 2 to a tr = 1-5, we find b tr ~ 24 for xh = 0.01 and b tT ~ 15 for xh = 0.1 M. 
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Finally, in Fig. 6 we present x*/xo as a function of the horizon radius Xh for b — 12.5, b — 21.125, b — 32 and 
b = 50. For b = 12.5 and b = 21.125 the monopole solutions show RN-type behaviour. This behaviour extends to 
the black hole solutions for b = 12.5 up to Xh < 0.17 and for b = 21.125 up to Xh < 0.07. The precise transition 
value must be extrapolated Jl(|. For larger values of Xh NA-type behaviour arises. For b = 32 and b = 50 regular 
and black holes solutions show NA-type behaviour. (For b = 105 and xt = 0.1 only the second minimum is left, 
as compared to b > 400 for the regular solutions Q.) For larger values of b and horizon radius Xh > 0.3 the 
analysis of black holes solutions needs special consideration H . 



V. CONCLUSION 



For gravitating monopoles and non-abelian black holes we have studied numerically the transition between 
the low and the high Higgs mass regime. Considering the approach of the solutions towards a limiting critical 
solution, the low Higgs mass regime shows RN-type behaviour, i.e. the critical solution where a horizon first 
appears corresponds to an extremal RN solution outside the horizon, while it is non-singular inside. In contrast, 
the high Higgs mass regime shows NA-type behaviour, i.e. the critical solution is an extremal black hole with non- 
abelian hair. In particular, analytical analysis suggests that RN-type behaviour should not occur for a < atr = 1.5 
corresponding to b > b tr ~ 25.6 Q. 

Our numerical analysis indicates two main results: 

• For asymptotically flat gravitating monopoles NA-type behaviour is still present at b = 30, considerably 
below the numerical value given by Lue and Weinberg 0] . We strongly expect that the NA-type behaviour 
persists until the analytical transition value is reached. But our numerical accuracy is insufficient to show 
this other than by extrapolation. 

• The analogous phenomenon occurs for black hole solutions. Here with increasing values of the horizon 
radius Xh, the transition occurs for decreasing values of b. 



VI. APPENDIX 



We employ a collocation method for boundary-value ordinary differential equations developed by Ascher, Chris- 
tiansen and Russell pjj] . The set of non-linear coupled differential equations eqs. ( jToj ) -(|l3|) is solved using the 
damped Newton method of quasi-linearization. At each iteration step a linearized problem is solved by using a 
spline collocation at Gaussian points. Since the Newton method works very well, when the initial approximate 
solution is close to the true solution, the gravitating monopole solutions for varying parameters a and b are 
obtained by continuation. 

The linearized problem is solved on a sequence of meshes until the required accuracy is reached. For a particular 
mesh Xi — x\ < X2 < ■■■ < ^tv+i = x , where Xi and x are the boundaries of the interval, and hi = Xi + \ — Xi, h = 
maxjx.j^jv/ii, a collocation solution v (x) = («i, «2, ■■■Vd) is determined. Each component v n (x) S C mn ~ 1 [x i ,x ], 
is a polynomial of degree smaller than k + m n , where m„ is the order of the n-th equation, and k is an integer 
bigger than the highest order of any of the differential equations. The collocation solution is required to satisfy 
the set of differential equations at the k Gauss-Legendre points in each subinterval as well as the set of boundary 
conditions. 

When approximating the true solution u n (x) by the collocation solution v n (x), an error estimate in each 
subinterval ) is obtained from the expression 

\\u®(x) - v®{x)\\m = c n Mn +mn) ^)\h k + m -- 1 + 0(h k+m »- l+1 ) (23) 



I = 0,....,m n - 1 , n = l,..., <2, (24) 

where c n j are known constants. Also, using eq. (^4|) a redistribution of the mesh points is performed to roughly 
equidistribute the error. With this adaptive mesh selection procedure, the equations are solved on a sequence of 
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meshes until the successful stopping criterion is reached, where the deviation of the collocation solution from the 
true solution is below a prescribed error tolerance pd| ] . 

For the gravitating monopolc solutions, we typically specified the error tolerance in the range Kr 4 -1(T 6 , but 
the absolute errors of the collocation solutions normally turned out to be far better, e.g. in the range 10 -7 — 10~ 12 . 
The number of mesh points used in these calculations was typically about 350, with nearly 70% of the mesh points 
in the critical region of the two minima of iV(x), if present. 

We calculated the solutions on the finite interval [0, x a ] with x ranging from 10 to 10 4 . The solutions are 
highly independent of the size of the interval, as seen in Fig. 7, where the metric function N(x) is shown for the 
parameters a = 1.43596, b — 36.125 and three value of x . Indeed, as specified by their boundary conditions 
eqs. (p7|)-(p^|), the gravitating monopole solutions are asymptotically flat. 

In Fig. 5, where the evolution of the two minima of the metric function N(x) is shown for b — 36.125, the 
corresponding global error estimates for the function N(x) are also shown. 

For the black hole solutions the set of boundary conditions eqs. (p0|)-(p2|) is employed, and the calculations are 
performed on an interval [xh, x ], where Xh represents the horizon of the black hole. The error tolerances and the 
errors are comparable to those of the gravitating monopole solutions. 
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Figure captions 

Fig. 1 

The metric function N(x) is presented in the NA-type regime for b = 36.125 and a = 1.43482, 1.43550, 1.43583 
and 1.43596. 
Fig. 2 

The critical value a cr — 1 is presented as a function of b in the NA-type regime for monopole solutions and for 
black hole solutions with horizon radii Xh = 0.01 and 0.1. 
Fig. 3 

The ratio x*/xq is presented as a function of b in the NA-type regime for monopole solutions and for black hole 
solutions with horizon radii Xh = 0.01 and 0.1. 
Fig. 4 

The value of the metric function N(x) at the outer minimum at the minimal value of a where the inner minimum 
appears is presented as a function of 6, for monopole solutions and for black hole solutions with horizon radii 
x h = 0.01 and 0.1. 
Fig. 5 

The values of the metric function N(x) at the outer minimum and at the inner minimum are presented as a 
function of a close to a cr for b = 36.125. Also shown is the global error estimate for N(x). 
Fig. 6 

The ratio x*/xq is presented for black hole solutions as a function of Xh in the NA-type regime for the values 
b = 12.5, 21.125, 32, and 50, corresponding to = 5, 6.5, 8, 10. 
Fig. 7 

The metric function N(x) is presented in the NA-type regime for b — 36.125 and a = 1.43596, as obtained on 
the interval [0, x a ] with x = 10, 1000, 10000. (The three curves cannot be distinguished graphically.) 
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